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Abstract 

We study the spin-orbit coupling in the whole Brillouin zone for GaAs using 
both the sp 3 s*d 5 and sp^s* nearest-neighbor tight-binding models. In the T-valley, 
the spin splitting obtained is in good agreement with experimental data. We then 
further explicitly present the coefficients of the spin splitting in GaAs L and X 
valleys. These results are important to the realization of spintronic device and the 
investigation of spin dynamics far away from equilibrium. 
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1 Introduction 



Spin-orbit coupling (SOC) is the key ingredient to the semiconductor spin- 
tronic devices [1,2]. Most of the proposed schemes of electrical generation, 
manipulation and detection of electron spin rely on it. Complete understand- 
ing of the SOC is therefore of great importance. In the bulk zinc-blende-type 
semiconductor such as GaAs, it is well known that near the center of Bril- 
louin zone the zero field splitting caused by the SOC depends cubically on the 
wave-vector k due to the bulk inversion asymmetry [3,4] or linearly due to the 
structure inversion asymmetry [5,6]. There are few investigations of the SOC 
for the states away from the band edge. The ab initio band structure calcula- 
tion [7], diagonalizing of truncated k • p Hamiltonian [7] or nearest-neighbor 
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tight-binding (TB) model including the SOC [8,9] have been performed to 
obtain the spin-orbit splitting outside of the Brillouin zone center. For the 
states near other high symmetry points such as the L and X points, one can 
get the form of the splitting from the symmetry property [10], whereas the 
actual coefficients need to be further calculated. Recently the splitting of the 
L- valley in bulk GaSb and GaSb/AlSb quantum wells were calculated using 
an sp 3 s*d 5 nearest-neighbor TB model including the SOC [8]. It is shown that 
the splitting in GaSb L-valley exceeds 10 meV, an order of magnitude larger 
than the typical value in the T-valley. For GaAs the corresponding data are 
still not available. 

The lack of quantitative information of the SOC outside the Brillouin center 
is not crucial to the development of spintronics at the present stage, since the 
electrons in most of the study locate at the bottom of the T-valley. However, 
in real situation, the devices usually work under high electric field which can 
drive the electrons to the states far away from T point, or even further to 
other valleys such as L- and/or X- valleys. Therefore for the realization of the 
spintronic devices, the SOC in the whole Brillouin zone is essential. In the 
previous works on high field spin transport in GaAs [11,12], the coefficient of 
spin splitting in high valleys are approximated by other material such as GaSb 
due to the lack of the corresponding data for GaAs. In this report we present 
the spin splitting of GaAs for the whole lowest conduction band. Especially, 
we calculate the the coefficients of spin splitting in L- and X- valleys. 



2 Calculation and Results 

Our calculations are performed in sp 3 s* and sp 3 s*d 5 nearest-neighbor TB mod- 
els with the SOC. This method has been proven to be an effective approach 
in band structure calculations [9,13,14,15,16,17,18]. The parameters we use 
are adopted from the published literate [9,13,14,15,16,17,18]. The original pa- 
rameter sets have been optimized to fit the experimental data, such as the 
band edge and the effective mass. With the comparison to the experimental 
data, the spin splitting near the T-point is calculated as the benchmark of the 
suitability of these parameters for calculating the spin splitting. 

In Fig. 1 we show the spin splitting around T point for different momentums. 
In the figure the parameter sets are chosen from Refs. [15] and [9] for sp 3 s* and 
sp 3 s*d 5 respectively. One can see from the figure that for small momentum the 
splittings calculated from sp 3 s* and sp 3 d 5 s* approaches both increase with the 
momentum. When the momentum becomes large, the splitting is no longer a 
monotonic. The splittings along different directions have different behaviors. 
For the states near the T-point, the splitting can be described by 7$"2(k) cr with 
cr being the Pauli matrices. In the coordinate system of x = [100], y = [010] 
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Fig. 1. (Color Online) (a) Spin splitting in T-valley for different momentum along 
T-K (solid curves) and T-W (dashed curves) directions; (b) The corresponding SOC 
coefficient 7 vs. momentum along T-K (solid curves) and T-W (dashed curves) 
directions; (c) 7 vs. angle of the momentum with different amplitudes. Solid curves: 
k = 0.01; Dashed curves: k = 0.05; Chain curves: k = 0.1; Dotted curves: k = 0.2 
(27r/a), respectively. Red/Blue curves are the results from the sp 3 s* d 5 / sp 3 s* model. 

and z = [001], 

7 0(k) = 7 (*,(*; - kl), ky(kl - kl), k Z (kl - fcj)). (1) 

This gives rise to a splitting that varies cubically in k, i.e., AE ~ 7/c 3 . 
The splitting along different directions are different. In the x-y plane, the 
splitting is AE = 7& 3 | sin(2#)| where 9 is the angle between the momen- 
tum and (100)-axis. The exact value of 7 is still in debate. Different ap- 
proaches give various values range from 8.5 to 34.5 eV-A 3 . Overviews of these 
results are nicely listed in Refs. [19,20]. Theoretical calculations based on 
different approaches give quite different values. There are two kinds of ex- 
periments that can measure 7 value. One is through the direct measurement 
of the splitting using Raman scattering. Experiments based on this approach 
show that 7 is about 23.5 eV-A 3 in wide GaAs quantum well [21]. In asym- 
metric GaAs/AlGaAs heterostructure/quantum well, 7 is about 16.5 or 11.0 
eV-A 3 [22,23]. The other kind of measurement is through spin relaxation time 
or magneto-conductance. This kind of measurement is indirect since it de- 
pends on how to qualitatively calculate the spin relaxation time or magneto- 
conductance. Earlier works of this kind estimate that 7 is about 20-30 eV-A 3 
[24,25,26]. Recent calculation based on fully microscopic approach reveals that 
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the experiments in two-dimensional (2D) system can be explained by using 
much smaller 7 value [27]. In our calculation, value of this coefficient is cal- 
culated as 7(k) = AE / (2\Q(k)\) . The results for different momentums are 
shown in Fig. 1(b) and (c). One can see from the figure that for k < 0.047r/a, 
7 is almost a constant that is independent on magnitude and the angle of mo- 
mentum. For the parameters we use, sp 3 s* and sp 3 s*d 5 give 7 = 17.0 and 23.9 
eV-A 3 near the T-point, respectively. Both are close to the experimental data 
from Raman scattering. This good agreement between the theoretical result 
and the experimental data shows that the parameters we use can be applied to 
study the spin splitting for whole Brillouin zone. One can see from the figure 
that both models predict that the value of 7 decreases with the increase of 
momentuni" 2 "! Moreover, angle dependence of 7 becomes remarkable for large 
momentum. Thus 7 is no longer a constant for large momentum. Practically 
the SOC described by Eq. (1) with constant 7 is a good approximation for the 
state not far away from equilibrium since even for high-carrier-density samples 
the value of 7 at Fermi surface is only a few percents smaller than the value 
at k — 0. However, when electrons are driven far away from the T point, Eq. 
(1) is expected to over-estimate the SOC. 

We now turn to the states in L- and X-valleys. The spin splittings for different 
momentums are plotted in Figs. 2 and 3(a) respectively. One can see from the 
figures that, for the same amount of momentum variation from the valley 
bottom, the spin splitting in L/X- valley is much larger than that in T- valley. 
For the states near the X-point, the SOC is in the form [10], (3 x fl x (k) ' a 
with 

(3 x n x (k)=(3 x (k x ,-k y ,0) . (2) 

Near the L- valley bottom, the SOC reads /^fi^k) • cr. In the new coordinate 
system spanned by [110] (x'-axis), [112] (y'-axis) and [111] (z'-axis) vectors, 
n L (k) reads [10] 

(3 L Q L (k)=(3 L (k' y ,-k' x ,0) . (3) 

In the above two equations, k represents the momentum vector measured 
from L/X- valley bottom. These couplings give the splitting linear to the first 
order of momentum around the valley bottoms. The corresponding splitting 
coefficients (3l/x — AE(k)/2k for different momentum are plotted in Figs. 
2(b) and 3(b) respectively. Similar to that of T- valley, these coefficients are 
constants near the valley bottoms. In the AT-valley, the values of f3 x obtained 
from sp 3 s*d 5 and sp 3 s* models are close to each other, i.e., f3 x = 0.059 and 
0.046 eV-A, respectively. However, in the L-valley, (3l determined from these 
two models are quite different. For sp 3 s*d 5 , @l = 0.26 eV-A; while for sp 3 s*, 
Pl = 0.047 eV-A. This profound difference implies that the d orbit plays an 
important role in the spin splitting in the L-valley of the lowest conduction 
band. This is because the symmetry imposes a cLorbital component in the 
L-valley and sp 3 s*d 5 model can account this symmetry more accurate. It has 

2 This explains the small value obtained in 2D system [27]. 
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Fig. 2. (Color Online) (a) Spin splitting in L-valley for different momentums along 
L-K (solid curves) and L-W (dashed curves) directions; (b) The corresponding SOC 
coefficients (3l vs. momentum along L-K (solid curves) and L-W (dashed curves) 
directions; (c) (3l vs. angle of the momentum with different amplitudes. Solid curves: 
k = 0.01; Dashed curves: k = 0.05; Chain curves: k = 0.1; Dotted curves: k = 0.2 
(2-7r/a), respectively. Red/Blue curves are the results from the sp 3 s*d 5 /sp 3 s* model. 



been revealed that the inclusion of the d orbit greatly improves the accuracy 
of the effective mass in L- valley [9,28,29,30,31]. Therefore, in our opinion the 
spin splitting determined by sp 3 s*d 5 model in L-valiey is also more reliable 
than that by sp 3 s* model. 



3 Conclusion 



In conclusion we study the SOC in the whole Brillouin zone for GaAs using 
both sp 3 s*d 5 and sp 3 s* nearest-neighbor TB models. For the parameter sets 
we use, the spin splittings calculated from both models are in good agreement 
with experimental data in the T-valley. We then further explicitly present the 
coefficients of the spin splitting in the L and X valleys. These results are useful 
for understanding the spin dynamics far away from equilibrium. 
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Fig. 3. (Color Online) (a) Spin splitting in X-valley for different momentums along 
X-K (solid curves) and X-W (dashed curves) directions; (b) The corresponding 
SOC coefficient fix vs. momentum along X-K (solid curves) and X-W (dashed 
curves) directions; (c) fix vs. angle of the momentum with different amplitudes. 
Solid curves: k = 0.01; Dashed curves: k = 0.05; Chain curves: k = 0.1; Dotted 
curves: k = 0.2 (2n/a), respectively. Red/Blue curves are the results from the 
sp 3 s* d 5 1 sp 3 s* model. 
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